Experimental verification of democratic particle motions by direct imaging of glassy 

colloidal systems 
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We analyze data from confocal microscopy experiments of a colloidal suspension to validate pre- 
dictions of rapid sporadic events responsible for structural relaxation in a glassy sample. The 
trajectories of several thousand colloidal particles are analyzed, confirming the existence of rapid 
sporadic events responsible for the structural relaxation of significant regions of the sample, and 
complementing prior observations of dynamical heterogeneity. The emergence of relatively compact 
clusters of mobility allows the dynamics to transition between the large periods of local confinement 
within its potential energy surface, in good agreement with the picture envisioned long ago by Adam 
and Gibbs and Goldstein. 

PACS numbers: 61.20.Ne (Structure of simple liquids), 82.70.Dd (Colloids) 



A complete understanding of the molecular underpin- 
nings of glassy relaxation (the dramatic dynamical slow- 
ing down that arises when a liquid is rapidly cooled be- 
low its melting point avoiding crystallization), remains 
a major challenge in condensed matter physics [IHB]. 
As long ago as 1965, Adam and Gibbs [7j proposed an 
appealing picture that accounted for this enormous in- 
crease in relaxation timescales within a narrow temper- 
ature window. They suggested that the dynamics of a 
glass-forming supercooled liquid proceeds by means of 
cooperatively rearranging regions (CRR) whose size and 
relaxation timescale grow considerably as temperature is 
lowered, giving the decrease in configurational entropy of 
the system [TH7] . This description suggests that a super- 
cooled liquid at low temperatures can be decomposed in 
independently relaxing compact subsystems (the CRR) 
whose molecules attempt to change configuration, but 
which can only undergo a transition when they rearrange 
in a concerted manner. Thus each of them (in the words 
of Adam and Gibbs 7J) surmounts, essentially simultane- 
ously, the individual barrier restricting its arrangement. 
In this picture, it is expected that each region of a su- 
percooled liquid should be practically "frozen" in a given 
portion of configuration space for large times (larger as 
temperature decreases given the growing size of the re- 
gions and thus of the number of molecules involved in 
the rearrangement) and then will relax (asynchronously 
and independently of other regions) by having a burst of 
mobility characterized by the sharp emergence of a com- 
pact cluster of mobile particles[71IH]. Hence, at any given 
time the system would present dynamics that would vary 
significantly from one region to another: the dynamics 
should then be heterogeneous in space [IHB]. 

The validity of such a heterogeneous scenario has 
been confirmed both experimentally and computation- 
ally, since the existence of dynamical heterogeneities [5]- 



HT] has been detected. Simulations of model glassy sys- 
tems have shown that the more mobile particles are 
not homogeneously distributed in space but arranged in 
(non-compact) clusters [TU]. The time scale for the mo- 
tion of these more mobile particles is t* , a timescale close 
to the structural cv-relaxation time, r Q ; r a is calculated as 
the time scale when the incoherent self intermediate scat- 
tering function has decayed to 1/e. These results have 
also received experimental support in colloidal suspen- 
sions (experimental models for glassy relaxation) [T21ll3| . 

More recently 18J, computational studies have deter- 
mined that within any dynamically heterogeneous region 
of the system, the relaxation is not gradual but also het- 
erogeneous in time, since the a relaxation is almost ex- 
clusively governed by rapid sporadic events character- 
ized by the emergence of relatively compact clusters of 
mobile particles (termed as "democratic" clusters or d- 
clusters [H]). These events trigger transitions between 
local metabasins (MB, basins of attraction of the poten- 
tial energy surface comprising a group of similar closely- 
related structures or local minima [3J IT8H2U] where the 
system is confined for long times). These cooperatively 
relaxing units or d-clusters have been identified in molec- 
ular dynamics simulations of different glassy systems like 
a binary Lennard- Jones [TS] system, supercooled water 
PTj and amorphous silica [22] and represent natural can- 
didates for the CRR proposed by Adam and Gibbs [7] . A 
recent inhomogeneous mode-coupling theory of dynam- 
ical heterogeneity has related them to the (fractal) ge- 
ometrical structures carrying the dynamical correlations 
at timescales commensurable with that of the a relax- 
ation (more compact than the open-like structures ex- 
pected at much shorter timescales) [23]. Additionally, a 
recent experimental and computational work in a glassy 
polymer provided indirect experimental support to the 
MB-MB transitions and d-clusters [21] • However, exper- 
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iments with molecular glasses lack the level of resolution 
necessary to directly observe them and thus, no direct ex- 
perimental information has verified the existence of such 
events up to date. 

In this work we study a colloidal suspension (an ex- 
cellent experimental model of glassy systems with parti- 
cles big enough to be directly observed by confocal mi- 
croscopy) and we shall demonstrate that detailed track- 
ing of particle motions is indeed able to detect the afore- 
mentioned kind of events. 

We analyze the data of Refs. [13l [14], taken from 
confocal microscopy experiments of colloidal samples. 
The colloids are sterically stabilized colloidal poly- 
methylmethacrylate) with diameter d = 2.36 /im and 
a polydispersity of ~ 5%. They are dyed with rho- 
damine and suspended in a density-matching and index- 
matching solvent mixture of cycloheptylbromide and de- 
calin. In this solvent, the colloidal particles possess a 
slight charge, and exhibit a glass transition at a volume 
fraction <j> s» 0.58. A confocal microscope rapidly ac- 
quires three-dimensional images once every 10 — 20 s. 
The images are post-processed to locate particle centers 
with an accuracy of 0.03 /im in x and y and 0.05 /im 
in z. Due to the difficulty of identifying particles near 
the edges of the images, the useful data are within a re- 
gion of size L x = 67 //m, L y — 62 /im, and L z = 9 /im, 
corresponding to several thousand particles. The volume 
fractions are determined by counting the particles within 
a subvolume, and are known to within ±0.01 with the 
uncertainty mainly due to the uncertainty of the particle 
diameter (±0.01 /im). For further experimental details, 
see Refs. [13 E]. 

Since at any given time a large system would consist of 
several different CRR, we divided the experimental sys- 
tem into 6 subsystems or portions £, each one with an 
increasing number of colloidal particles N(£) (see the in- 
set in Fig. [lja)). All portions have the same depth L z 
and center (L x /2,L y /2, L z /2). Portion £ comprises the 
particles that were initially (t = 0) within the bound- 
aries of the corresponding rectangular prism of length 
Ac ' (£/6) and height L y ■ (£/6). For 4> — 0.56, the num- 
ber of particles within each portion is: = 1) = 77, 
AT(2) = 310, A(3) = 703, AT(4) = 1255, N(5) = 1962 
and N(6) = 2759 particles. We shall present results for 
(f) = 0.56 unless otherwise indicated, but similar results 
were obtained for <j> = 0.53, <fi = 0.52 and = 0.46. 

To identify MBs we employed the following "distance 
matrix" (A 2 ) function [25] : 



A 2 (t',0 = ^El r ^')-r^")| 2 , (1) 

where r.j(i) is the position of particle i at time t. A 2 (t' , t") 
gives the system normalized squared displacement in the 
time interval (t' , t"). 
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FIG. 1. (a) Distance matrix A 2 (i', t") for portion £ = 1. The 
gray level corresponds to values of A 2 (t', t") that are given to 
the right of the figure. Units are /im 2 . Inset: for the analysis, 
we followed the colloidal particles that were initially (t — 0) 
within the boundaries of rectangular prisms: £ = 6 (blue, the 
complete system), magenta, green, yellow, orange and £ = 1 
(red, smallest region) respectively. All £ have the same depth 
L z . Colloidal particles are not shown, (b) Averaged squared 
displacement 5 2 (t, 9) for all £. Each series (analysis over each 
£) is shifted by 0.5 /im 2 respect the former one. For compar- 
ison we included the corresponding average values of 5 2 (t, 6) 
over all times, (S 2 (t,8)) (dashed colored lines, also shifted). 
The value of 6 is 72 s. Inset: 7 = j<5f(5706 s, 6)-(5 2 (t, 9)} s \/a s 
vs £, where a is the standard deviation in S 2 and 5706 s is the 
time of the largest average squared displacement. Subscript £ 
means that the function was evaluated for molecules in £. A 
maximum for £ = 2 is observed, (c) The function m(t, 6) for 
£ = 1 and its average value over all times (dashed black line) . 

A plot of A 2 as a function of t' and t" can be seen in 
Fig. [lja) for <fi = 0.56. These results are typical for all 
studied 4> (as an example, in Fig. [2] we show an equiv- 
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FIG. 2. Similar to Fig. [lfa), but for (ft = 0.46 and for a 
(rectangular prism) region of approximately the size of £ = 1. 
At this 0, iV(l) = 44 and t* = 300 s. 

alent plot for (ft = 0.46). The darker the shading, the 
smaller the distance between the configurations at times 
t! and t". From this figure we can learn that the dy- 
namics of this portion is quite heterogeneous in time in 
that it stays for a significant time relatively close to one 
region in configuration space, dark square-like regions, 
before it finds a pathway to a new region. The value of 
A 2 (t',t") within a MB is around 0.02 /zm 2 as compared 
to values much larger than 0.04 /im 2 if the system is in 
different MBs (see legend on the right of the figure). If 
there were no MBs, the plot would show a dark shadow 
at the diagonal t" — t' and a gradual decrease in shading 
perpendicular to it, as it would be seen at low (ft (simi- 
larly to the case of structural glasses at high temperature 
[TB] ) and/or large systems; compare Fig.JTJa) (eft = 0.56) 
to Fig. H ((ft = 0.46). These fi gures demonstrate that 
the system spends large amounts of time exploring the 
local MB, and only occasionally moves on to a neighbor- 
ing MB. Indeed, we can see that this trajectory resides 
within a MB for times much larger than t* w 1000 s, the 
maximum in the non-Gaussian parameter a.<i (t) [T3l I15j . 
We also point out that from Fig. [lja) it is evident that 
the time for a MB-MB transition is quite short, on the or- 
der of 70 s, which thus corresponds to about 7% of t* . In 
the lower volume fraction data of Fig. [2j the transitions 
are also rapid although slightly less distinct. 

In Fig. W[b) we also show, for all £ and same time 
interval, 5 (t,9), the averaged squared displacement of 
the particles within a time interval 6 (solid curves). This 
function is defined as 

8 2 {t,9) = A 2 (t,t + 9) (2) 
1 N 

= ^Ei r ^)- r *( < + ^ 2 ■ ( 3 ) 
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FIG. 3. Distribution functions for £ = 2 and 9 = 72 s. The 
black curve is 4nr 2 G B (r, 6), the self-part of the van Hove func- 
tion, and the gold curve is a Gaussian with the same value of 
{r 2 (0)). The crossing point of these two curves at r ss 0.23 fim 
is used as a threshold to identify the democratically moving 
particles. The blue and magenta curves are the average of 
4nr 2 G a (r, t,t + 6) for different values of t in which the system 
is inside a MB (see text for details), and a transition (TR) 
from t = 5706 s to t + 6. 

A comparison of 5 2 with the distance matrix shows 
that S 2 is showing pronounced peaks exactly when the 
system leaves a MB. Thus we see that changing the MB is 
indeed associated with a rapid significant particle motion 
as measured by S 2 . Also included are the average values 
of 8 2 {t,9) for all £ over all times, (5 2 (t,9)), represented 
by dashed lines. It is clear that the larger systems have 
a lower relation between fluctuations in S 2 at the MB 
transitions and its average value: This ratio in Fig. [ljb) 
for the peak at time t = 5706 s is maximum at £ = 2 
(see inset), thus providing an indication of the size of the 
MB-MB transition event. 

To understand the motion of the particles when the 
system leaves a MB we have calculated the function 
4:Trr 2 G s (r,t,t + 9), the distribution of displacement r 
of the particles for a given time difference 6 — 72 s. 
(Note that the average of 4irr 2 G s (r,t,t + 6) over t gives 
47rr 2 G s (r, 9), the self-part of the van Hove function). An 
average of this distribution is shown in Fig. [3] (blue curve) 
for different values t within a MB [t/s = 378, 2898, 3474, 
4626, 6210 and 7506 for the case of Fig. [lja)]. Also in- 
cluded is the self-part of the van Hove function (black 
curve) and we can see that both curves are (within the 
noise of the data) identical, thus showing that in a MB 
the system moves basically the same as on average. We 
also show the distribution (magenta curve) for t = 5706 s 
in which the system is about to leave a MB. For this value 
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FIG. 4. Actual size and position of democratic particles 
for: (left) a MB-MB transition from t = 5706 s to t + 6 and 
(right) t = 7290 s to t + 9 inside a MB. The data are for £ = 2, 
= 72 s, $ = 0.56. 

of t the distribution is clearly displaced to the right with 
respect to 4irr 2 G s (r, 8), showing that in this time regime 
the motion of the system is much faster than on average. 
Thus we can conclude that the peaks of the 5 2 are not 
due to the presence of a few fast moving particles, but 
instead to a "democratic" movement of many particles. 

To further demonstrate that this is indeed the case 
and to explore the spatial distribution of mobility, we 
have defined as "democratic" all those particles that 
in the time interval 8 = 72 s have moved more than 
fth = 0.23 /im (a value very close to 0.25 /xm, the size 
of the cage formed by particles surrounding a single par- 
ticle [2]), and denote the fraction of such particles by 
m(t,8). We take the value of r t h as the second inter- 
section between knr 2 G s (r, 8) and a Gaussian (see Fig. |4j 
gold curve) with the same value of mean squared dis- 
placement, MSD = (r 2 (8)) (however, other threshold 
choices yield similar results). In Fig. [ljc) we have in- 
cluded the fraction m of democratic particles for £ = 1 
as a function of time (vertical bars) . A comparison of this 
data with 8 2 shows that m is indeed large whenever S 2 
increases rapidly. This fraction is on the order of 20 % of 
the particles and thus significantly larger than one would 
expect from 4irG s (r,8) if one integrates this distribution 
from rth to infinity and which gives 0.051. In turn, Fig. [3] 
(left) shows the 3D location of the democratic particles 
involved in a typical MB-MB transition for £ = 2 (we 
show the case for the event at t = 5706 s in Fig. [T] but 
other cases display similar results). We can see that the 
particles are not homogeneously distributed in space (as 
is indeed the case in Fig. [1] (right) for t — 7290 s inside a 
MB) but arranged in a relatively compact cluster. 

In summary, the experimental results shown in this 
work for a colloidal system provide direct validation to 
the picture of glassy relaxation previously shown by MD 
simulations of several glassy systems [THl HU [221 123] : 
The dynamics spends large times confined within a 
metabasin, interspersed with rapid bursts in mobility 
characterized by the emergence of relatively compact 
clusters of democratic particles which trigger the struc- 



tural or a relaxation. Systems closer to the glass transi- 
tion have more distinct metabasin transitions that occur 
more infrequently (compare Fig. [lja) to Fig. [2]). Thus, 
this behavior conforms to the scenario put forth long time 
ago by Adam and GibbsjT] and Goldstein [5]. 
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